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ABSTRACT 

Based on updated nonlinear, convective pulsation models computed for several val- 
ues of stellar mass, luminosity and metallicity, theoretical constraints on the AT— band 
Period-Luminosity (PLk) relation of RR Lyrae stars are presented. We show that for 
each given metal content the predicted PLx is marginally dependent on uncertainties 
of the stellar mass and/or luminosity. Then, by considering the RR Lyrae masses sug- 
gested by evolutionary computations for the various metallicities, we obtain that the 
predicted infrared magnitude Mk over the range 0.0001< Z <0.02 is given by the rela- 
tion M K =O.568-2.O711ogP+O.O871ogZ-O.7781ogL/L , with a rms scatter (7^=0.032 
mag. Therefore, by allowing the luminosities of RR Lyrae stars to vary within the range 
covered by current evolutionary predictions for metal-deficient (0.0001 < Z <0.006) 
horizontal branch models, we eventually find that the infrared Pcriod-Luminosity- 
Metallicity {PLZ K ) relation is M K =0.139-2. 071(logP+0.30)+0.1671ogZ, with a to- 
tal intrinsic dispersion of <jk =0.037 mag. As a consequence, the use of such a PLZk 
relation should constrain within ±0.04 mag the infrared distance modulus of field and 
cluster RR Lyrae variables, provided that accurate observations and reliable estimates 
of the metal content are available. Moreover, we show that the combination of K and 
V measurements can supply independent information on the average luminosity of 
RR Lyrae stars, thus yielding tight constraints on the input physics of stellar evolu- 
tion computations. Finally, for globular clusters with a sizable sample of first overtone 
(RR C ) variables, the reddening can be estimated by using the PLZk relation together 
with the predicted My— logP relation at the blue edge of the instability strip (Caputo 
et al. 2000). 

Key words: globular clusters: distances - stars: evolution - stars: horizontal branch 
- stars: oscillations - stars: variables: RR Lyrae 



1 INTRODUCTION 

Together with Classical Cepheids, RR Lyrae stars are pop- 
ular standard candles for estimating stellar distances in the 
Milky Way and to Local Group galaxies. Cepheids are the 
typical Population I standard candles, while RR Lyraes are 
generally observed in globular clusters (GCs), thus provid- 
ing a widely adopted method to calibrate the luminosity of 
the main-sequence turn-off and to constrain the age of these 
very old stellar systems. 

In the last few years distances based on the RR Lyrae 
properties uncorked a flourishing literature. However in spite 
of these efforts, we are still facing the problem that dis- 
tance determinations based on different empirical methods 
present discrepancies that are larger than the estimated in- 
trinsic errors (see e.g. Carretta et al. 2000), thus suggest- 



ing the presence of poorly known systematic errors. On the 
other hand, theoretical predictions concerning the horizontal 
branch (HB) luminosity at the RR Lyrae gap are affected by 
uncertainties in the input physics of stellar models that do 
not allow us to nail down accurate estimates of the absolute 
magnitude Mv(RR) of RR Lyrae stars (Cassisi et al. 1998; 
Caputo 1998; Castellani & Degl'Innocenti 1999; Caputo et 
al. 2000). 

The pros and cons of the various empirical methods cur- 
rently adopted to determine RR Lyrae distances both in the 
field and in GCs have been widely discussed in the litera- 
ture and they are not further addressed here. However, it is 
already known that near infrared (NIR) observations, and 
in particular K— band photometry, can overcome several 
deceptive errors affecting the RR Lyrae distance scale. As 
a matter of fact, K magnitudes present some indisputable 
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advantages when compared with usual B or V magnitudes: 
i) a much smaller dependence on interstellar extinction, ii) 
a smaller dependence on stellar metallicity, and finally iii) 
smaller pulsational amplitudes which allow us to reach a 
superior precision on time averaged magnitudes with a rel- 
atively small number of observations. 

The NIR observational scenario for RR Lyrae stars has 
been further enriched by the empirical evidence brought out 
by Longmore, Fernley & Jameson (1986) and by Longmore 
et al. (1990, hereinafter L90) who found that cluster RR 
Lyraes do obey a rather precise PLk relation. On the ba- 
sis of RR Lyrae K photometry in eight Galactic GCs, L90 
found that the intrinsic scatter of individual PLk relations 
was dominated by observational errors, as well as that all the 
relations have a slope cLMk / 'dlogP ~ —2.2, within the em- 
pirical uncertainties. The occurrence of a PL relation much 
steeper in the K than in the V ~ band is obviously due to 
the fact that K magnitudes present a stronger dependence 
on the effective temperature than the optical ones (see also 
Nemec, Nemec & Lutz 1994). As far as the intrinsic accu- 
racy of the PLk relation is concerned, it might be taken as 
evidence that the scatter of stellar masses and/or luminosi- 
ties among RR Lyrae stars marginally affects the correlation 
between K magnitudes and periods (L90). 

To calibrate the empirical PLk relations, L90 provided 
three zero-points according to the absolute distances of field 
RR Lyrae stars, as measured with different Baade-Wesselink 
methods. In all cases, a mild dependence on metallicity 
was suggested, with the [Fe/H] coefficient ranging from 
0.04 to 0.08 mag dex~ x . However, L90 also found that cal- 
ibrations with the IR flux method (Fernley et al. 1989; 
Skillen et al. 1989) give RR Lyrae luminosities ~0.14 mag 
brighter than predicted by the IR Baade-Wesselink method 
(Jones, Carney, & Latham 1988, hereinafter JCL88; Liu & 
Janes 1990, hereinafter LJ90). More recently, Jones et al. 
(1992, hereinafter J92) presented a detailed analysis of 18 
field RR Lyrae stars with distances estimated according to 
the IR surface brightness version of the Baade-Wesselink 
method and with metallicities ranging from solar down to 
[Fe/H]=— 2.2. The authors concluded that the PLk relation 
appears largely insensitive to the stellar metal abundance. 

Summarizing, the most recent empirical metal-free 
PLk relations available in the literature are those provided 
by J92 

M K = -0.88(±0.06) - 2.33(±0.20) logP 

and by Skillen et al. (1993) 

M K = -1.07(±0.10) - 2.95(±0.10) logP, 

while, if the metallicity term is taken into account, the em- 
pirical period- luminosity-metallicity (PLZk) relations are 

M K = -0.647 - 1.72 log P + 0.04 [Fe/H], 
M K = -0.76 - 2.257 log P + 0.08 [Fe/H], 

and 

M K = -0.72(±0.11)-2.03(±0.27) log P+0.06(±0.04) [Pe/if] 

by L90, LJ90 and J92, respectively. We add that the abso- 
lute calibration suggested by L90 according to the IR flux 
method distances of two field RR Lyrae stars yields quite 
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Figure 1. Top panel: predicted My magnitudes for fundamen- 
tal (filled symbols) and first overtone (open symbols) pulsators 
at fixed mass (0.65Mq) and chemical composition (Y=0.24, 
Z =0.001), but with three different luminosities (see labeled val- 
ues). The models cover the entire pulsation region, from the blue 
to the red edge of the instability strip and the periods of FO pul- 
sators are fundamentalised. The solid lines display the mean vi- 
sual magnitude of the pulsators at fixed luminosity. Bottom panel: 
predicted Mk magnitudes for the same pulsators. The solid lines 
show the PLk relations at fixed luminosity. 

brighter magnitudes, namely 

Mif.-o.so = -0.24 + 0.06 [Fe/H], 

where M K ,-o.so is the absolute infrared magnitude at 
logP=-0.30. 

At variance with such a large body of observational 
work, we still lack a detailed theoretical study of the physi- 
cal scenario underlying the PLk relation. In this paper, we 
wish to investigate this topic on the basis of an extensive 
grid of nonlinear, convective pulsation models constructed 
by adopting chemical compositions, luminosities and stellar 
masses typical of field and cluster RR Lyrae variables. In 
Section 2 we present the adopted sets of models and discuss 
the dependence of the predicted if— magnitudes on period 
and metallicity. The theoretical PLZk relation is discussed 
in the same section, together with the comparison with cur- 
rent empirical data for RR Lyrae stars in M3. Final remarks 
and future developments close the paper. 



2 PREDICTED PLZ K RELATION 

As quoted before, empirical studies of the PLk relation 
suggest that both the evolutionary (i.e., luminosity) effects 
and the spread of stellar masses inside the instability strip 
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Figure 2. Same as in Fig. 1, but with fixed luminosity 
(logL/Lo=1.61) and two different mass values, namely O.625M0 
(triangles) and 0.675Mq (circles). 



marginally affect that relation. To assess the plausibility of 
such a scenario, we computed three sequences of pulsating 
models for fixed mass (0.65Mq) and chemical composition 
(Y=0.24, Z=0.00T), but with three different luminosities 
(logL/L0=1.51, 1.61, 1.72) and effective temperatures rang- 
ing from the blue to the red edge of the instability region. 
The theoretical framework adopted to construct the non- 
linear convective models has been presented in a series of 
papers (Bono & Stellingwerf 1994; Bono, Caputo, & Mar- 
coni 1995; Bono et al. 1997b; Caputo et al. 2000) and it is 
not further discussed here. The predicted luminosities have 
been transformed into the observational plane by adopting 
the bolometric corrections and the color-temperature rela- 
tions provided by Castelli, Gratton, & Kurucz (1997a, b). 

Figure 1 shows the predicted My (top panel) and Mk 
(bottom panel) magnitudes for fundamental (F) and first 
overtone (FO) pulsators, as a function of the pulsation pe- 
riod. Note that the FO models have been fundamental- 
ized by adding 0.127 to their logarithmic period, while the 
predicted magnitudes are magnitude-weighted means over 
the pulsation cycle. It is worth noting that mean magni- 
tudes can be derived by adopting two different methods, 
i.e. magnitude- weighted (m,j) or intensity- weighted < rrij >. 
However, the NIR amplitudes are relatively small and there- 
fore the discrepancy between synthetic (K) and < K > 
magnitudes is less than 0.01 mag for both F and FO pul- 
sators, in close agreement with empirical studies (see e.g. 
L90 and Fernley 1993). Concerning the visual band, the 
difference (V)— < V >, as well as the difference between 



statky and mean magnitudes, depend on the symmetry of 
the light curve, generally increasing as the amplitude of pul- 
sation increases. As shown by Bono, Caputo & Stellingwerf 
(1995), the differences for the first-overtone mode are always 
smaller than for fundamental mode, and among F pulsators 
the discrepancies are significantly increasing moving from 
the almost sinusoidal light curves near the red edge to the 
saw-tooth ones near the blue edge of the instability region. 
Also this finding agrees with empirical evidence (Fernley 
1993). 

A glance at the data plotted in Fig. 1 brings out the 
mechanism of the K band "degeneracy". For any fixed lu- 
minosity level, the pulsation period increases from the blue 
to the red edge of the instability region and the My magni- 
tude remains, as expected, almost constant, except for some 
faster fundamental pulsators appearing fainter according to 
the light curve effect already discussed in Bono et al. (1995). 
At any fixed luminosity the dependence of My on period is 
dMy/dlogP « —0.14, whereas the dependence of Mk is 
steeper, and indeed dMx/dlogP w —2.1. This effect is due 
to fact that when moving from the blue to the red edge of the 
instability strip the period increases and the pulsators be- 
come brighter, due to the strong dependence of the K-band 
bolometric correction on effective temperature. An increase 
in the luminosity, for any fixed effective temperature, causes 
an increase in the period according to dlogP/dlogL ~0.84, 
while the visual magnitude closely follows the expected shift 
in luminosity as dMy/dlogL ~ —2.5. On the other hand, 
a change in the luminosity level causes a variation in both 
K magnitudes and periods in such a way that the PLk re- 
lation of bright pulsators fairly match the relation for the 
fainter ones (see bottom panel of Fig. 1). 

Figure 2 shows My (top panel) and Mk (bottom panel) 
magnitudes for models with the quoted chemical composi- 
tion, fixed luminosity (logL/Z/0 = 1.61), but with two differ- 
ent mass values (0.625 and 0.675Mq) taken as a reasonable 
estimate of the expected mass dispersion among RR Lyrae 
stars with logZ ~ —3.0. We find that the PLk relation is 
only marginally affected by the spread in the stellar mass, 
and indeed a variation of 8 percent introduces, at fixed pe- 
riod, a variation in Mk of ~0.06 mag. 

Based on the results presented in Fig. 1 and Fig. 2, 
we find that for a given chemical composition the effect on 
the PLk of a dispersion in the stellar mass of ±4% is al- 
most negligible (~ ±0.03 mag), and the effect caused by a 
dispersion in luminosity is significantly smaller than in the 
visual relation. Even if we account for a dispersion as large as 
A log L = ±0.1, the PLk relation would give the absolute K 
magnitude (i.e. the distance modulus) with an accuracy bet- 
ter than ±0.1 mag, which is substantially smaller than the 
analogous dispersion in the visual band, i.e. AMy ~ ±0.25 
mag. 

Theoretical predictions on the PLk relation do de- 
pend on the adopted bolometric corrections and color- 
temperature relations. Recent detailed investigations sug- 
gest that predictions based on current atmosphere models 
are in good agreement with empirical data for both dwarf 
and giant stars (Bessell, Castelli, & Plez 1998). The same 
outcome does not apply to RR Lyrae stars, and indeed dis- 



The magnitude such a star would have were it not pulsating. 
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Table 1. Input parameters for the pulsation models. 



z 


Y 


M/Mq 


log L/Lq 


0.0001 


0.24 


0.75 


1.72, 1.81 


0.0004 


0.24 


0.70 


1.61, 1.72, 1.81 


0.001 


0.24 


0.65 


1.51, 1.61, 1.72 


0.006 


0.255 


0.58 


1.55, 1.65 


0.01 


0.255 


0.58 


1.51, 1.57 


0.02 


0.28 


0.53 


1.41, 1.51, 1.61 



crepancies of the order of 100-200 K have been found be- 
tween temperature estimates based on optical and NIR col- 
ors (Clementini et al. 1995; Cacciari et al. 2000). Fortunately 
enough, the V-K colors seem only marginally affected by 
systematic uncertainties, and indeed Clementini et al. found 
that predicted temperatures are, at fixed V-K, on average 50 
K cooler than observed ones. To estimate the impact of this 
uncertainty on Mk magnitudes we assumed an error of ±50 
K on the temperature and we found that predicted V-K col- 
ors for both metal-poor and metal-rich structures change by 
±0.03 mag. As a consequence, the uncertainty introduced by 
the (V-K)-temperature calibration on Mk seems negligible 
when compared with evolutionary effects. 

To investigate in more detail the NIR period-luminosity 
relation, theoretical evaluations have been extended to the 
metallicity range 0.0001 < Z < 0.02, by adopting for each 
given value of Z a helium content fixed according to a Galac- 
tic helium-to-metal enrichment ratio AY/AZ ~2.5. For each 
chemical composition we adopted the mass of the Zero Age 
Horizontal Branch (ZAHB) model (Bono et al. 1997a; Cas- 
sisi et al. 1998) located in the middle of the instability strip 
(logT e ~3.85), while several luminosity levels have been ex- 
plored both to cover the different theoretical ZAHB lumi- 
nosities predicted in the recent literature (Caputo et al. 
2000, and references therein) and to account for the expected 
evolutionary effects. Table 1 summarizes the input param- 
eters of the computed pulsating models, while Table 2 and 
Table 3 give the details of fundamental and first overtone 
pulsators, respectively^. 

The top panel of Fig. 3 shows the predicted Mk mag- 
nitudes for the entire set of models plotted as a function of 
period. We find that Mk is tightly correlated with period, 
metallicity and luminosity according to the relation 

M K = 0.568-2.071 log P±0.087 log Z-0.778 log L/L Q (1) 

with a rms scatter ctk=0.016 mag. 

The accuracy of eq. (1) in reproducing the predicted 
Mk magnitudes is disclosed in the bottom panel of Fig. 
3, where the magnitudes are scaled to Z=0.001 and 
logL/Lo=1.61, according to eq. (1). Taking into considera- 
tion the uncertainty of ~0.03 mag due to a mass dispersion 
by ±4% with respect to the values in Table 1, we even- 
tually find that the total intrinsic dispersion of eq. (1) is 
± 0.032 mag. As a result, we conclude that over the en- 
tire explored metallicity range a variation of AlogL=±0.1 
moves the predicted K magnitudes by ±0.08 mag, while a a 
change in metallicity of 0.3 dex causes a change of AK ~0.03 

t Table 2 and Table 3 list magnitude-weighted quantities over the 
pulsation cycle. Intensity-weighted magnitudes can be requested 
to marcella@na.astro.it. 
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Figure 3. Top panel: predicted Mk magnitudes of F (filled cir- 
cles) and FO pulsators (open circles) as a function of period (data 
in Table 2 and Table 3, respectively). Bottom panel: same as the 
top panel, but with the Mk magnitudes scaled to the same metal- 
licity (Z=0.001) and luminosity (logL/L =1.61). The solid line 
shows the relation predicted by eq. (1). In both panels the periods 
of FO pulsators are fundamentalised. 



mag. This evidence implies that RR Lyrae infrared magni- 
tudes are excellent standard candles to derive distances of 
field and cluster RR Lyrae with accurate estimates of metal 
abundances. Furthermore, as we will discuss in the follow- 
ing, the milder sensitivity of the PLk relation to the stellar 
luminosity supplies the opportunity for testing theoretical 
predictions on the ZAHB luminosity. 

As an example, Fig. 4 shows the observed K— logP 
and V— logP diagrams of RR Lyrae stars in the Galactic 
globular cluster M3 (data collected by L90 and by Car- 
retta et al. 1998, respectively). The slope of the predicted 
PLk relation (dashed line) appears in reasonable agreement 
with observations, even though there is a significant disper- 
sion in the empirical data which seems to be dominated 
by observational errors, as suggested by L90. However, new 
and more accurate data are necessary to assess whether 
this dispersion is intrinsic or caused by observational er- 
rors. For the sake of the discussion, let us take into account 
the mean of all the visual magnitudes, i.e. V =15.65±0.10 
mag, and the averaged infrared magnitude at logP=— 0.30, 
i.e. K =14.65±0.05 mag. By assuming for M3 an iron-to- 
hydrogen ratio [Fe/H]=-1.34 (Carretta & Gratton 1997), 
i.e. \ogZ=— 3.04, the data in Table 2 and Table 3 (see also 
the top panel of Fig. 1) provide the predicted averaged Mv 
magnitude of the pulsators, as a function of the adopted lu- 
minosity. At the same time, under the same assumptions on 
metallicity and luminosity, eq. (1) yields the predicted Mk 
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Table 2. Periods and absolute mean magnitudes (magnitudc-wcighto 
and effective temperatures. 
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magnitude at logP=— 0.30. The absolute visual and infrared 
magnitudes are listed in Table 4 [column (2) and column (3), 
respectively], together with the corresponding visual (fiv) 
and infrared (hk ) distance moduli of the RR Lyrae stars in 
M3 [column (4) and column (5), respectively]. 

Inspection of the data in Table 4 shows that an increase 
in the luminosity of AlogL=0.21 causes an increase in the 
distance modulus of 0.52 mag in the visual but only 0.16 mag 
in the infrared. Since the reddening towards M3 is negligible 



(Dutra & Bica 2000) and can be set to zero, the average lu- 
minosity < log L/Lq >rr of RR Lyrae stars can be found 
as the luminosity for which the distance moduli in V and 
K assume the same value (fio = fix = /iv). Thus, interpo- 
lation in Table 4 leads to < logi/L Q > fli{ = 1.68±0.05 for 
a mean distance modulus /ix = 15.03±0.07 mag, where the 
error on the distance modulus also includes the observed dis- 
persion of infrared data. We wish to notice that this value 
of the M3 distance modulus is in close agreement with that 
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Table 3. Same as Table 2, but for first overtone pulsators. 
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(fiv = 15.00±0.07 mag) derived by Caputo et al. (2000) 
on the basis of the comparison between the observed dis- 
tribution in the V— logP plane of M3 c— type variables and 
the predicted blue edge for first overtone pulsators, as given 
by the same pulsating models adopted in the present paper. 
This agreement confirms the internal consistency of our pul- 
sational scenario, as well as that the reddening towards M3 
is almost zero. 

In this context it is worth emphasizing that the sensitiv- 
ity of < log L/Lq >rr to errors in the reddening determina- 
tion has a rather negligible effect for the true distance mod- 
ulus derived by infrared observations. As a matter of fact, if 
we assume for M3 a reddening E(B — V) ~0.05 mag. With 
A V /E(B - V)=3.1 and A K /E(V - *0=0.13 (see Cardelli, 
Clayton & Mathis 1989), the apparent visual and infrared 
moduli in Table 4 should be decreased by 0.16 mag and 
0.02 mag, respectively, and the condition fio,v = Mo, if leads 
now to a higher luminosity (< log L/Lq >_h_r=1.76±0.05), 
but for a quite similar true distance modulus, i.e. Ho,k = 
15.07±0.07 mag. On a general point of view the depen- 
dence of < log L/Lq >rr on the adopted reddening is 
8 log L/8E(B — V) ~1.6 dex mag -1 , whereas the sensitivity 
of the true distance modulus based on infrared observations 
is only 5li ,k/5E(B - V) ~0.9. 

Finally, we wish to notice that our distance evaluation 
is in close agreement with the M3 distance estimated by L90, 



Table 4. Predicted M v and M K (at logP=-0.30) magnitudes 
for pulsators with logZ=— 3.04, as a function of the adopted lu- 
minosity (see text). The visual (fJ-v) an d infrared ((j-k) distance 
moduli of RR Lyrae stars in M3 are derived according to the ob- 
served averaged magnitudes V=15.65 mag and K(— 0.30)=14.65 
mag. The last column gives the difference between visual and 
infrared distance moduli. 



log L/Lq My 

1.51 1.045 

1.61 0.797 

1.72 0.523 



-M"if,-0.30 

-0.250 
-0.328 
-0.413 



14.605 
14.853 
15.127 



14.900 -0.295 
14.978 -0.125 
15.063 +0.064 



i.e. jj,K = 15.00 ± 0.04 ± 0.15, where the uncertainties refer 
to errors in the zero-point and in the slope, respectively. As 
a consequence, our theoretical approach to the PLZk rela- 
tion seems to support the absolute calibration based on the 
IR flux method. A detailed comparison between individual 
cluster distances based on empirical methods and theoreti- 
cal predictions, as well as a thorough analysis of the error 
budget, will be addressed in a forthcoming paper. 



3 DISCUSSION 

In the previous section it has been shown that, if one re- 
lies on theoretical pulsation results, the predicted PLZk 
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Figure 4. Top panel: V— logP diagram of RR Lyrae stars in 
M3 from data collected by Carretta et al. (1998). The solid line 
shows the average visual magnitude (15.65 mag). Bottom panel: 
K— logP diagram of RR Lyrae stars in M3 from data collected 
by L90. The dashed line has the slope predicted by eq. (1), while 
the solid line is the average infrared magnitude (14.65 mag) at 
log P=— 0.30. In both panels the periods of c— type variables are 
fundamcntaliscd. 



Figure 5. Top panel: luminosity (solar units) versus mctallic- 
ity of ZAHB models located in the middle of the RR Lyrae gap, 
as given by Bono et al. (1997a, filled circles) and Cassisi ct al. 
(1998,1999, open circles). The asterisk marks the average lumi- 
nosity of RR Lyrae stars in M3, as derived in the present pa- 
per. Bottom panel: predicted K magnitudes at logP=— 0.30 for 
the two sets of models, according to cq. (1). The solid line is 
eq. (2), while the dashed lines refer to the PLZjf relations with 
\ogZ < -3.0 and logZ > -3.0. 



relation of RR Lyrae stars can supply accurate distance 
estimates to globular clusters with reliable infrared obser- 
vations. In fact, this relation appears only marginally de- 
pendent on the exact value of the luminosity as well as on 
current uncertainties on the iron-to-hydrogen content. More- 
over, by combining K and V magnitudes one can derive in- 
dependent constraints on the average bolometric magnitude 
of RR Lyrae variables. 

In this context, we show in Fig. 5 (top panel) the 
logL(ZAHB)— logZ relations based on model computations 
by Bono et al. (1997a, filled circles) and by Cassisi et al. 
(1998, 1999, open circles), which cover current uncertainties 
on the ZAHB luminosity, together with the average luminos- 
ity of M3 RR Lyrae stars (asterisk), as derived in the previ- 
ous section. We find that the ZAHB luminosities predicted 
by Bono et al. (1997a) agree quite well with the expected 
average overluminosity (SlogL ~0.04) of actual RR Lyrae 
stars when compared with the ZAHB luminosity. On the 
other hand, the HB models computed by Cassisi et al. (1998, 
1999) using the most updated physics do predict ZAHB lu- 
minosities that are systematically too bright. This seems 
to strengthen the evidence based on independent estimates 
(Bono et al. 2000; Caputo ct al. 2000) that the most updated 
physics could be still affected by systematic uncertainties. 

However, the bottom panel of Fig. 5 shows that the 
current uncertainty on ZAHB luminosity has a rather weak 



effect on K magnitudes at fixed period. Then, by adopting 
the two logL(ZAHB)— logZ calibrations plotted in Fig. 5 as 
a measure of the luminosity dispersion among actual vari- 
ables with a given metallicity, we find that the predicted K 
magnitude of RR Lyrae stars with metal content in the range 
0.0001< Z <0.006 is correlated with period and metallicity 
as 

M K = 0.139 - 2.071(logP + 0.30) + 0.167 log Z, (2) 

with a total intrinsic dispersion of ctk=0.037 mag, including 
the above uncertainty on luminosities and mass variations 
by ±4% on the values listed in Table 1. If we also wish 
to account for the nonlinearity of the logL(ZAHB)— \ogZ 
relationship and for the change of the slope at \ogZ ~ —3.0, 
then the least-squares solutions for models with logZ < —3.0 
and logZ > — 3.0 supply 

M K = 0.058 - 2.071(log P + 0.30) + 0.141 log Z (3) 

and 

M K = 0.213 - 2.071(log P + 0.30) + 0.196 log Z (4) 

respectively, with <tk=0.035 mag. 

It is worth mentioning that for globular clusters with a 
sizable sample of c— type variables, the visual distance mod- 
ulus may be derived by matching the observed RR C distri- 



© 0000 RAS, MNRAS 000, 000-000 



8 G. Bono et al. 



bution in the V-\ogP plane with the predicted first overtone 
blue edge (FOBE). According to Caputo et al. (2000), the 
Period-Luminosity-Metallicity relation for FOBE pulsators 
is 

My = -0.1 78 - 2.255 log P + O.f 5f log Z, (5) 

with a total intrinsic dispersion of ov=0.065 mag. It follows 
that the comparison of K and V magnitudes with theo- 
retical predictions gives the opportunity to derive reliable 
estimates of the cluster reddening. Fortunately enough, the 
coefficients of the metallicity term in eq. (2) and eq. (5) are 
approximately the same, and therefore they cancel out any 
metallicity effect on RR Lyrae colors. 

Of course, the results presented in this paper depend 
on the reliability of the theoretical scenario concerning both 
pulsational models and bolometric corrections. To quanti- 
tatively assess their intrinsic accuracy is obviously a dif- 
ficult task, owing to the rather sophisticated physics in- 
volved in the modeling of stellar pulsation and stellar at- 
mosphere. However, at least for the pulsation framework, 
we have shown that two quite different observables based on 
the same pulsating models, namely the PLZk relation given 
in the present paper and the Mv-logP relation for FOBE 
pulsators provided by Caputo et al. (2000), consistently pro- 
vide similar values for the M3 distance modulus. Moreover, 
in a recent paper it has been shown that the same pulsation 
models can closely reproduce the details of the observed light 
curve of RR Lyrae stars (Bono, Castellani & Marconi 2000) . 
As a whole, these results can be taken as further evidence 
that the theoretical framework we are constructing in the 
last years is approaching the behavior of actual pulsating 
stars. 

Finally, let us note that even though K observations 
present several significant advantages when compared with 
optical bands, it is limited how deep one can reach from the 
ground due to the high level of thermal background emission. 
De facto NIR data are only available for Galactic RR Lyrae 
stars, and even with the largest ground based telescopes, the 
observations can only be extended out to a few Local Group 
galaxies. However, NGST could change this picture substan- 
tially with its superior NIR capability, and bring most Local 
Group galaxies within reach. The PLZk relation does pro- 
vide a method to determine very precise distances to these 
galaxies based on well understood physics. Therefore, it will 
supply a very important independent check on the Cepheid 
distance scale, as well as reliable calibrations of secondary 
distance indicators. 
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